Boundary conditions and defect lines in the Abelian sandpile model 
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We add a defect line of dissipation, or crack, to the Abelian sandpile model. We find that the defect 
line renormalizes to separate the two-dimensional plane into two half planes with open boundary 
conditions. We also show that varying the amount of dissipation at a boundary of the Abelian 
sandpile model does not affect the universality class of the boundary condition. We demonstrate 
that a universal coefficient associated with height probabilities near the defect can be used to classify 
boundary conditions. 

PACS numbers: 05.65.+b,45.70.-n 
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I. INTRODUCTION 

The Abelian sandpile model (ASM) was introduced by 
Bak, Tang, and Wiesenfeld as a model of self-organized 
criticality yj. This well-known model was designed to 
demonstrate how simple rules can drive a system to a 
critical point, and thus produce power laws, without any 
fine-tuning of parameters. It has thus been used to ex- 
plain power laws in a wide range of systems — see 0,0 for 
a review. Since the ASM was first introduced, a number 
of variations on the model have been introduced — see Q 
for a review. However, the original ASM still provides a 
simple, important, and robust model for the generation 
of power laws. 

The ASM is defined on a square lattice. Each site a 
of the lattice has a height variable, h a , which can be any 
integer from 1 to 4, inclusive, where h a represents the 
number of grains of sand at that site. At each timestep, 
a grain of sand is added to a random site of the lattice. 
After the addition of the grain, any site with more than 
four grains of sand is unstable, and collapses, losing four 
grains of sand, while each of its four neighbors gains one 
grain of sand. Unstable sites are repeatedly collapsed, 
until every site is stable — i.e. no site has more than four 
grains. Then, the next timestep, another grain is added, 
and the entire process is repeated. 

The original ASM is spatially homogenous (except for 
the boundaries, which break translational invariance), 
and most modifications of the sandpile model have kept 
this feature. However, here we consider the effects of a 
crack, represented by a defect line, along which grains of 
sand can be lost; in other words, along which the num- 
ber of grains is not conserved. In previous studies, dis- 
sipation was added to the bulk of the ASM (not break- 
ing translational invariance), and was shown to take the 
ASM off its critical point jE 0, ■ Our defect line of 
dissipation breaks translational invariance, and we show 
that it causes the two-dimensional plane of the ASM to 
renormalize into two half planes with open boundary con- 
ditions. This shows that cracks in the ASM are highly 
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relevant, and essentially cleave the sandpile into separate 
pieces. We demonstrate this by looking at the universal 
coefficient associated with the modification of unit height 
probabilities at large distances from the defect, and at 
the correlation function between unit height variables on 
opposite sides of the defect. The Green function for an 
ASM with a defect line is calculated in section ITTT1 and 
results for the height probabilities and correlations are 
presented in section llVl 

For most models of interest in condensed matter 
physics, the bulk properties can be studied with the 
boundary playing little to no role. For example, the two- 
dimensional Ising Model is often studied on a torus, so as 
to eliminate boundary effects. However, this is not possi- 
ble for the ASM. In the bulk of the ASM, the number of 
grains of sand is conserved during each toppling. If this 
was true for all sites, then eventually we would reach a 
state where topplings continued without end. The ASM 
thus needs sites with dissipation — that is, sites where the 
number of grains is not conserved. The most natural way 
to do this is with open boundary conditions; sites at the 
open boundaries become unstable when they have more 
then four grains (just as in the bulk), but have only three 
neighbors to send grains to, and send the fourth grain 
"off the edge," removing it from the system. Since this 
dissipation is necessary for a well-defined sandpile, the 
boundary plays a crucial role in the ASM, even when we 
are focused at points in the bulk. Correlation functions 
far from the boundary are independent of the boundary 
conditions, just as in other condensed matter statisti- 
cal mechanical models; but the presence of dissipation 
somewhere in the ASM (e.g. at the open boundary) is 
necessary for the model to be well-defined. 

We consider the effects of varying amounts of dissi- 
pation along a boundary, and show that any amount of 
dissipation at the edge results in the open boundary uni- 
versality class. The Green function is calculated in sec- 
tion and results for the height probabilities and cor- 
relations are presented in section IVTI 

These results are intuitively reasonable, since dissipa- 
tion should be relevant in regions of the ASM where the 
particles have no other way to leave the ASM. However, 
it was also possible that such modifications could have re- 
sulted in new, as yet undiscovered, boundary conditions 
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or defect states. For example, Bariev, and McCoy and 
Perk, added a line defect of modified bond strengths to 
the Ising Model, and found that they were able to contin- 
uously vary the dimension of the spin operator along the 
defect by varying the defect bond strength [g, |9( . This 
continual variation occurred despite the fact that the 
Ising Model only has three conformally invariant bound- 
ary conditions. 

The ASM has been associated with a conformal field 
theory (CFT). While CFT's are generally well under- 
stood, the ASM is a logarithmic CFT (LCFT) (specifi- 
cally, the c = — 2 CFT) , many aspects of which are still 
not well understood [T(|- In particular, our understand- 
ing of the boundary states of LCFT's is still fragmentary, 
and recent results on boundaries of the c = — 2 CFT have 
been partially contradictory [Til IT2I Il3l IT3. Il5j . Connec- 
tions between the LCFT bounary states, and the ASM 
boundary states were made in but the ASM rep- 
resentation of the some of the c = — 2 LCFT boundary 
states is still unknown. Modifications to the ASM such 
as those described in this paper, and searches for other 
boundary conditions, could help eludicate these relation- 
ships. Our results provide some evidence that the open 
and closed boundary conditions are the only possibilities 
for the ASM, although it is still possible that further cal- 
culations in this vein could uncover new boundary con- 
ditions. 

The identification of boundary states as closed, or 
open, or in some new, as yet undiscovered class, uses ar- 
guments from CFT that the coefficients associated with 
the falloff of expectation values (height probabilities) at 
large distances from the defect should be universal. Our 
results both use this expected universality, and confirm 
it, since we find, for example, that the coefficient is unaf- 
fected by varying a free parameter corresponding to the 
amount of dissipation at the boundary. This confirma- 
tion, while expected, is valuable, given the anomalous 
and unsettled nature of boundary LCFT associated with 
the ASM. This is a particularly important point in light 
of recent arguments that use this universality to argue 
that the four height variables in the ASM must corre- 
spond to different fields in the corresponding CFT [l7|. 



II. THE FORMALISM 

Dhar pointed out ASM is highly analytically tractable 
because of its Abelian nature — the same state results 
whether grains of sand are added first at site a and then 
at site 6, or first at site b and then at site a 18]. This 
is the basis of a well-established formalism for analyzing 
the ASM — see 01 f° r a review. We only give a quick 
coverage of the essential points here. 

It is useful to first generalize the above description of 
the ASM, to allow for more complicated topplings. The 
dynamics of the model are described by a toppling ma- 
trix, A a b, where a and b label sites of the lattice. The 
dimension of A is equal to the number of sites in the lat- 



tice, so A becomes infinite-dimensional as the size of the 
lattice goes to infinity. We say that site a is unstable if 
its height h a is greater than A aa . If site a is unstable, 
then every height changes by hb — > hb — A a b (including at 
the site b = a). We have the standard ASM, with open 
boundary conditions, if A a b is 4 when a = 6, -1 when a 
and b are nearest neighbors, and otherwise. 

Dhar showed that the states of any sandpile, given cer- 
tain general conditions on the form of A, are divided into 
transient states, which occur with probability zero after 
long amounts of time, and recurrent states, which all oc- 
cur with equal probability after long amounts of time. 
The number of possible recurrent configurations is given 
by det(A) Eg. 

Furthermore, Majumdar and Dhar also showed how 
to find the probability for a site to have height one, 
and the joint probability for two sites to both have 
height one (as well as other, more complicated proba- 
bilities) [20I |2]| . The toppling matrix is modified by 
removing specific bonds, and changing the toppling con- 
dition at certain sites. For example, if we want to force 
site a to have height 1 , we change the toppling matrix so 
that A aa — 1, and remove three of the bonds to neigh- 
boring sites (setting A a & = A{, a = for those bonds). 
With this modified toppling matrix, A', site a is now 
guaranteed to have height 1, and det(A') gives the num- 
ber of recurrent configurations with h a = 1. While A 
and A' are infinite-dimensional matrices (for an inifnite 
lattice), B = A' — A is outside of a 4x4 submatrix. So 
det(A')/det(A)=det(l + BA~ r ) is an easily computable 
4 by 4 matrix determinant, which gives the probability 
that, in a randomly chosen recurrent configuration, the 
site a will have height 1. The same process, with a dif- 
ferent (8 by 8) matrix B, can be used to find two-point 
correlations of height 1 variables. 

This process requires us to calculate the Green function 
G = A -1 . The Green function has long been known for 
the standard ASM, where A is simply the lattice Lapla- 
cian |23. However, in the following sections we will be 
dealing with different toppling conditions, and so will 
need to calculate the Green function for these new A's. 



III. GREEN FUNCTION FOR THE DEFECT 
LINE 

We introduce a defect line (or crack) through the mid- 
dle of the ASM, allowing dissipation to take place along 
the defect, and not just along the open boundary condi- 
tions. We take the lattice to be size M x [2L— 1), with the 
x-dimension taking on the values i = 0, 1, ... , (M — 1), 
and the y-dimension taking on the values j = —{L — 
1), —(L - 2), . . . , (L - 2), (L - 1). We take open bound- 
ary conditions along all edges, and put the defect along 
the line j = 0. Along this line, the height variable can 
take on the values 1, 2, . . . , (4 + m), where m > 0. A 
site along the defect topples if its height is greater than 
(4 + m). When it topples, it sends one grain to each of 
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its four neighbors, and to grains of sand are dissipated 
(i.e. disappear from the sandpile). 

When to is a positive integer, the theory has its most 
obvious physical interpretation, but the theory can be 
modified to give a sensible interpretation for any rational, 
positive, value of m 0- If in each toppling, ci grains are 
toppled, and C2 grains sent to each neighbor, where C\ 
and C2 are integers, then the ratio of grains dissipated to 
grains moved, to/4 «-> (ci/(4c2)) — 1, can be any rational 
integer. 

The toppling matrix A is the same as for the standard 
ASM, except that A aQ = 4 + m for sites a along the 
defect. When to = it becomes the standard ASM. The 
toppling matrix can be written as 



A ( V 



A (2 > 

33' 



= d^A^+S^Alp .with 
2 if i = %' 

— 1 if % = i' ± 1 , and 
otherwise 

2 ifj=j'^0 

to + 2 if j = j' = 

-1 ifj=j'±l 

otherwise 



(1) 
(2) 

(3) 



Since A is Hermitian, if we find all of its normalized 
eigenvectors, we can easily invert it. Suppose that the 
eigenvectors of A are ep t g, with eigenvalues \p. p and x 
are two-dimensional vectors and the number of possible 
values of p is equal to the dimension of A, which is in 
turn equal to the number of sites in the lattice. Then 



A- 1 - 



(4) 
(5) 



The form of A in equation2]implies that the eigenvectors 
of A factorize into eigenvectors of A^ and . 

We thus want the eigenvectors of A^ 2 ). (The eigenvec- 
tors of A^ 1 ' are not only simpler, but immediately follow 
from the eigenvectors of A^ 2 ) , by setting to = 0.) j and j' 
range from — (L — 1) to (L — 1), so A*- 2 ) has 2L — 1 eigen- 
vectors. The eigenvectors fall in three classes. There are 
(L — 1) oscillatory eigenvectors that are antisymmetric 
about j — 0, and have momenta p evenly spaced be- 
tween and 7r (p — rnr/L, n 6 2, 1 < n < (L - 1)). 
There are another {L — 1) oscillatory eigenvectors that 
are symmetric about j = 0, and have momenta p in the 
range < p < 7r, where the p solve a transcendental 
equation; in the limit L — > oo these momenta p also be- 
come equally spaced between and n. Finally, there is 
one exponentially decaying eigenvector, symmetric about 
3=0. 

Since A is Hermitian, we can immediately obtain its 
inverse from these eigenvectors. The sums over the two 
oscillatory sets of eigenvectors each produce integrals in 
the limit L — ► oo, M — ► oo, using the Euler-MacLaurin 
formula. The last, exponentially decaying, eigenvector 
produces a single, discrete contribution to the Green 
function. Writing the Green function as a sum of the 
contributions from the three classes of eigenvectors gives 



G(i,j,i',f) = J2G^(i,j,i',f) (6) 

a=l 

G^(i,j,i',f) = l ~G {i-i',j-j')-\Go{i-i',j+3') (7) 
G^(i,j,i',f) = l -G t) {i-i l 1 \ 3 \-\ 3 l \) + 8G^{i-i l 1 \ 3 \ + \ 3 l \) (8) 

-cm+ij'i) 



where we have defined 



K= 2 + +i . (10) 



Gn is the well-known bulk Green function 221: 



g a i)= [ 27! — f 27r dPl cos ^ gggCgjjj z 1 fn) 

J 2tt J 2tt 4 — 2 cos — 2 cospj 1 ' 



We have also defined 
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6G^(i,j) 



m 
~~2 



) G^ (i, j) + m (G^ (i, j-l)- G^ (i, j + 1)) 
I(g( 2q )( 1 ,j-2) + G( 2q )(* ! j + 2)) +Cm 



2 

2-rr j__ r2w 



' j7r dp., cos(pii) cos(pjj) 



dpi 
2tt 



27r 4 2 cos pj — 2 cos p., to 2 + 4 sin 2 p^ 
dpj 1 



sin 2 pj 



2ir 2 - cos pi - cospj m 2 + 4 sin 2 p 3 



(12) 
(13) 

(14) 



We want the behavior of G(i,j,i,j') for | j \ + | j' | 
large. The expansion of the bulk Green function Gq(0, j) 
for large j is well-known : 



r 



using the methods outlined in section [H] We find that 
the probability for a site a distance j from the defect to 
have height 1 is 



Go(0 ! i)^-il gO-)-i(2 + ^log2 



247rj 5 



(15) 

where 7 = 0.577... is the Euler-Mascheroni constant. 
We also need the behavior of <5G( 2a ' (0, j) for j large. 
The integral over pj in equation ED can be done exactly, 
and making the substitution z — e lpi gives a contour 
integral around the unit circle. The integrand has two 
poles inside the unit circle, but these give contributions 
which either decay exponentially with j, or are indepen- 
dent of j, neither of which affects our height calculations; 
so these contributions can be dropped. The algebraic j- 
dependence comes from the branch cut in the integrand, 
running from z = 3 — \/8 to z = 3 + v°j which gives 



G (2a) (0,j) 
/(*) 



1 r z 3 — 1 

-V \ dz -/(*) 

"K J3-V8 z - 1 
1 



(16) 



V-z 2 + 6z-l \ {z 2 -l) 



(17) 



where V indicates that we take the principal part of 
the integral. We can use this to find the behavior of 
C 2a '(0, j) for large j, by separating out the contribu- 
tions from z near 1, and expanding in a Laurent series in 
j. This then gives the expansion of <5G( 2 '(0, j): 



5G ^(0,j) = ^lo g j + ^ + ^^ + ..., (18) 

47T TO7TJ 4o7rm z j z 

dropping terms independent of j. 



IV. HEIGHT PROBABILITIES FOR THE 
DEFECT LINE 

Now that we have the Green function for the defect 
line, we can use it to calculate the height correlations, 



prob (h 



1 



2(tt - 2) 



1 



1 

4f 



1 



(id) -) ^3 ^ 1 4^2 2mj 3 

(19) 

The constant term, 2(w — 2) / w 3 , is the bulk probability 
for a site to have height one, first calculated in [2(J] . It was 
also found in [2(j that the correlation function between 
two height one operators is 



prob (fyi^) = h (i2th 
2(tt-2) x 2 



= 1 



1 - 



2r 4 



r = - 12) 2 + (h - J2) 2 



, where (20) 
(21) 



Thus, the height one operator is a dimension 2 oper- 
ator. Based on the identification of the ASM as a con- 
formal field theory 0, |2^|, the coefficient of 1/j 2 , in 
expectation values of dimension 2 operators a distance 
j from a boundary, is expected to be a universal num- 
ber characteristic of the boundary condition [23 |. And, 
in fact, this coefficient of +1/4 in equation 1191 is exactly 
the coefficient seen for the height one probability at large 
distances from an open boundary condition, as shown by 
Brankov, Ivashkevich, and Priezzhev |2fij . This indicates 
that, upon renormalization, the defect line becomes an 
open boundary. 

It is only sensible to talk about conventional boundary 
conditions at the defect line, if the two half-planes on ei- 
ther side of the defect have somehow been separated. Ev- 
idence that the defect renormalizes to separate the half- 
planes can be seen by looking at correlation functions of 
points on opposite sides of the defect. If, upon renor- 
malizing the defect, the two sides of the defect were still 
"connected," we would expect that height variables on 
opposite sides would still fall off as 1/r 4 , since the height 
one operator has dimension 2. (Calculations of corre- 
lation functions along boundaries by Ivashkevich have 
shown that the height one operator also has dimension 2 
along open boundaries HiJ.) However, we find that 
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prob (h(ij) = 1, h(i-j) = l) - prob = l) prob = l) = 

= ( ( tt3 } ) (~8^ + °(f)) 



While the height variable is a dimension two operator, 
its correlations across the defect fall off as 1/r 6 . The 
coefficient of the f/r 4 term renormalizes to zero, and the 
1/r 6 term is non-universal, depending continuously on 
in. We thus conclude that the defect renormalizes to 
generate two separate half-planes with open boundary 
conditions. 

This is physically reasonable. Adding dissipation 
throughout the bulk of the ASM is known to be rele- 
vant, driving the system off criticality |EIE0- More re- 
cently, adding dissipation in the bulk was identified with 
adding the integral of a dimension variable (the loga- 
rithmic partner of the identity) throughout the bulk |23j . 
ft would thus appear that the local addition of dissipa- 
tion should be represented by a dimension zero operator, 
which would mean that adding dissipation along a de- 
fect line should be relevant, as we have found here. In 
short, cracks in the ASM cleave the plane into discon- 
nected regions with open boundary conditions. Similar 
separation with relevant perturbations along a defect oc- 
curs in other models — see, for example |27j . However, 
these results were not inevitable; as already noted, line 
defects added by Bariev, and McCoy and Perk, to the 
fsing Model, resulted in a continual range of defect lines, 
along which the dimension of the spin operator could be 
continuously varied, despite the fact that the Ising Model 
only has three possible boundary conditions. 



GREEN FUNCTION FOR THE MODIFIED 
BOUNDARY CONDITIONS 



i = 0, 1, . . . , (M — 1), and the y-dimension taking on the 
values j = 0, 1, . . . , (L— 1). We impose a modified bound- 
ary condition along j = 0, and open boundary conditions 
on the other three edges. (In the end, we take the lim- 
its L — > oo and M — > oo, so our results should anyway 
be insensitive to the boundary conditions on these three 
edges.) 

We allow the height variable on the boundary, j = 0, to 
take on values I, 2, . . . , 6. Sites on the boundary become 
unstable when their height is greater than b, at which 
point they topple, giving one grain to each of their three 
neighbors, and dropping 6 — 3 grains off the edge. For 
6 = 3 this is the closed boundary condition, and for 6 = 4 
this is the open boundary condition. For 6 < 3, sand 
is generated with each toppling, rather than dissipated, 
creating the possibility of never-ending cycles of toppling. 
We therefore only consider 6 > 3. The system can be 
given a sensible interpretation for any rational value of 
6 > 3 5]. The toppling matrix between sites and 
(i',f) can be written 



where A' 1 ) was defined in equation[21 and 



(i) 



(23) 



2 if j = j' ^ 

(3) _ J b-2 if j=j' = 

ii' ~ 1 -1 ifj=j'±l 

otherwise 



(24) 



The identification of a line of dissipation with open 
boundary conditions brings up the question of whether 
any other universality classes of boundary conditions 
with dissipation are even possible. The open and closed 
boundary conditions are the most natural to impose on 
the ASM, but other boundary conditions than these two 
conventional ones can be written down. We create new 
boundary conditions by varying the amount of dissipa- 
tion along the boundary, and show that, regardless of 
the amount of dissipation, we stay in the open boundary 
universality class (so long as the amount of dissipation is 
nonzero — that is, so long as the boundary is not closed). 

It is convenient to change the dimensions of the lattice 
from those of section IIIII We take the lattice to be of 
size M x L, with the x-dimension taking on the values 



As with the defect, if we can find all the eigenvectors 
of A, we can easily invert it. A' 3 ', being L-dimensional, 
has L eigenvectors. When 3 < 6 < 5, A^ 3 ) has L eigen- 
vectors that are oscillatory functions of j, with momenta 
p, which satisfy a transcendental equation. In the limit 
L — ► oo, the momenta are evenly spaced over the range 
< p < it. When 6 > 5, A^ 3 ^ only has (L - 1) such 
oscillatory eigenvectors, and one last eigenvector that is 
exponentially decaying in j. 

In the Green function, the summation over oscillatory 
eigenvectors can be turned into an integral in the limit 
L — » oo, M — ► oo, with the Euler-MacLaurin formula. 
For 6 > 5, the single, exponentially decaying eigenvector 
produces a separate, discrete contribution to the Green 
function. The Green function is then given by 
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G(i,j,i',f) = G(i,j,i , ,j , )+e(b-5)G^(i,j,i , ,j') , where (25) 

d( . . , ;f) _ f 27T dm [*"d Pj coaQuCt-Q) 1 

{hJll ' J) J Q 2tt J 2tt 2-cosj)i-cosj)j (6 - 4) 2 + 1 + 2(6 - 4) cosft 

x [(6 - 4) sinfej) + sinfoO" + 1))] [(& - 4) sm( Pj f) + sm( Pl (f + 1))] , and (26) 

G°*>(i,j,f,f) = (6 ~J^_ 5) (4 - b)-i-i-\k . (27) 



We have defined 



k 

6{x) 



l + W^~) ' and 

1 if x > 
if x < 



(28) 
(29) 



At first sight, this equation for the Green function 



seems to indicate that G has a slope discontinuity at 
6 = 5. However, this is not the case. G is not smooth 
at b = 5, and expanding G as a function of b near 6 = 5 
shows that the combination G + 9(b — 5)G exp is actually 
smooth to all powers of (6—5). Inspection shows that G 
is also a smooth function of 6 for all other 6 in the range 
3 < 6 < oo (including 6 = 4). 

We need the Green function in two limits. First, for 
i = i' and j + f large, it is useful to write 



G(i, j, %', f) = G (i- i' t j -j')- G Q (i — + f) + 5G(i - i', j + f) 
SG(i,j) = 



2tt 



dpi 
2ir 



dpi 



cos(pji) 



an j 2tt 2 — cos j»j — cospj (6 — 4) 2 + 1 + 2(6 — 4) cos Pj 
x [sin P j cos(pjj) + (6 — 4 + cos^ ) sin(pjj)] . 



(30) 



(31) 



In (5G(0,j), we can do the integral over pi exactly, and 
then set z — e lpj . As before, the main contribution comes 
from the branch cut between z = 3 — y8 and 2=1. 
Expanding the integral near 2 = 1 gives 



SG(0,j) 



6 2 - 86 + 19 



7r(6-3)j tt(6-3) 2 j 2 27r(6-3) 3 j 3 



(32) 

We also need the expansion of the Green function along 
the defect — that is, for j = j' = and | i — i' |^> 0. Using 
similar methods as before, we find 



G(x 



,3=3 



0) 
6 2 



186+57 



7r(6-3) 2 a; 2 2vr(6 - A) A x A 



O (33) 



VI. HEIGHT PROBABILITIES FOR MODIFIED 
BOUNDARY CONDITIONS 

Using the Green function for modified boundary con- 
ditions, we can calculate unit height probabilities with 
the methods described in section ITTl We find that the 



probability for a site a distance j from the boundary to 
have height 1 is 



prob (h 



(id) 



2(tt-2) 



1 

4f 



2(6 -3) f 



(34) 

As discussed earlier, the coefficient of the 1/j 2 term is 
expected to be a universal number characteristic of the 
boundary condition [24(, and is equal to +1/4 for the 
open boundary condition [2^|. We see here that the co- 
efficient is +1/4, and independent of 6 for 6 > 3. This 
both confirms the expectation that the coefficient should 
be universal, and indicates that the boundary is in the 
open boundary universality class for any amount of dis- 
sipation (6 > 3). 

Note that the coefficient of 1/j 3 is non-universal, and 
diverges as 6 — > 3, indicating that 6 = 3 is a special point 
as we vary 6. 6 = 3 corresponds to the closed boundary 
condition, and it is already known that the coefficient of 
the 1/j 2 term is different (-1/4) for the closed boundary 
conditions; this is appropriate, since the closed and open 
boundary conditions are clearly in different universality 
classes |'25j. 

Boundary correlations along the j — boundary can 
be calculated, and contain no surprises. The correlation 



function between sites (i, 0) and («', 0) falls off as 1/ | i — i' | for all values of b 



prob {h {ifi) 
fl 



1, fyi',o) = l) - prob (/i (ii0 ) = 
(1 - c + c 2 )(-l + 3c - 4ci + c 2 ) 



(&-3)tt 



1) prob (/i(i/,o) 

2 



1 







(35) 



coefficient 




b 



4 5 6 7 8 9 10 



FIG. 1: Coefficient of the 1/r 4 term in the two-point correla- 
tion function 

We have defined c x = G(x —\ i — i 1 \,j=j' = 0). Equa- 
tions 1251271 can be used to find analytic expressions for 
c x , for x — 0, 1, 2. However, the expressions are long and 
not particularly enlightening, so are not presented here. 
The (absolute value of the) coefficient of the 1/ | i — i' | 4 
term is plotted in figure ^ It falls on " smoothly with 
increasing b. 

The coefficient of 1/ | i — i' | 4 in equation 1351 diverges 
as b approaches 3, reflecting the fact that b = 3 is a 
fixed point of the renormalization group flows, leading 
to non-smooth behavior in physical properties. However, 



the Green function, and height correlations calculated 
from it, are perfectly smooth as we vary b through 4. ft 
would appear that the RG flows take us from b = 3 to 
b = oo, and that b = 4 is not a fixed point of the RG 
flows. However, b = oo is in a sense the same as b = 4, 
in that both equally well represent the open boundary 
condition; if b = oo, the sites j — can hold an infinite 
number of grains, and never topple — the sandpile thus 
acts as if j = 1 was the boundary, with an open boundary 
condition, where grains fall "off the edge" to j = 0. 

We have shown that the addition of dissipation along a 
defect line separates the ASM into two half-planes, each 
with open boundary conditions. This brought up the 
question of whether there are other universality classes 
of boundary conditions, with varying amounts of dissipa- 
tion along the boundary. We find that any amount of dis- 
sipation along a boundary results in the open boundary 
universality class at large distances. Classes of bound- 
ary conditions were identified by the universal coefficient 
of the unit height probability, far from the boundary or 
defect. 
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